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ABSTRACT 

We present deep iKTnear-IR and Spitzer mid-IR observations of a large sample of spectroscopically- 
confirmed galaxies at z > 6. The sample consists of 51 Lya emitters (LAEs) at z ~ 5.7, 6.5, and 7.0, 
and 16 Lyman-break galaxies (LBGs) at 5.9 < z < 6.5. The near-IR images were mostly obtained 
with WFC3 in the F125W and F160W bands, and the mid-IR images were obtained with IRAC in 
the 3.6/im and 4.5/im bands. Our galaxies also have deep optical imaging data from Subaru Suprime- 
Cam. We utilize the multi-band data and secure redshifts to derive their rest-frame UV properties. 
These galaxies have steep UV continuum slopes roughly between /3 ~ — 1.5 and -3.5, with an average 
value of /? ~ —2.3, slightly steeper than the slopes of LBGs in previous studies. The slope shows 
little dependence on UV continuum luminosity except for a few of the brightest galaxies. We find 
a statistically significant excess of galaxies with slopes around j3 ~ — 3, suggesting the existence of 
very young stellar populations with extremely low metallicity and dust content. Our galaxies have 
moderately strong rest-frame Lya equivalent width (EW) in a range of -10 to -200 A. The star- 
formation rates are also moderate, from a few to a few tens solar masses per year. The LAEs and 
LBGs in this sample share many common properties, implying that LAEs represent a subset of LBGs 
with strong Lya emission. Finally, the comparison of the UV luminosity functions between LAEs and 
LBGs suggests that there exists a substantial population of faint galaxies with weak Lya emission 
(EW < 20 A) that could be the dominant contribution to the total ionizing flux at z > 6. 
Subject headings: cosmology: observations — galaxies: evolution — galaxies: high-redshift 



1. INTRODUCTION 

The epoch of cosmic reionization marks one of the ma- 
jor phase transitions of the Universe, during which the 
neutral intergalactic medium (IGM) was ionized by the 
first astrophysical objects. Measurements of CMB polar- 
ization ha ve shown that reioni zation began earlier than 
z ~ 10.6 (|Komatsu et al.ll2011[ ). Meanwhile, studies of 
Gunn-Peterson troughs in high-redshift quasar spectra 
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have located the end of reionization at z ~ 6 (|Fan et al.l 
2006). Therefore, objects at z > 6 are natural tools 
to probe this epoch. With recent advances of instru- 
mentation on Hubble Space Telescope (HST) and large 
ground-based telescopes, galaxies at z > 6 are now be- 
ing routinely found. Direct observations of the earliest 
galaxy format ion and the history of reionization a re fi- 
nally in sight ([Robertson et al .1120 lot iFinlatorl 120121 ). 

The first z > 6 galaxies were discovered to be Lya 
emitters (LAEs) at z ~ 6.5 using the narrow-band 
(or Lya) technique (|!LeLaLll2QQ2|; iKodaira et all 120031 : 
iRhoads et al.l 120041 ). This technique has been an ef- 
ficient way to find high-redshift galaxies due to its 
high success rate of spectroscopic confirmation. Three 
dark atmospheric windows with little OH sky emission 
in the optical are often used to detect distant galax- 
ies at z ~ 5.7, 6.5, and 7. So far more than 200 
LAEs have bee n spectroscopically co n firmed at these 
redshifts (e.g., iTaniguchi et al.l [20051: llv e ^jap 120061: 
Shi masaku et al. 2006; Kashikawa e t al.l I2006L 120111 : 
Hu et all 120101 : lOuchi et al.l 120101 : IRhoads et"aLl 



The narrow-band technique is also be ing used to search 
for higher redshift LAEs at z > 7 (e., 
Tilvi et allfeOltt IClement et aLlBQl 



Hibon et al. 2010; 



2012 



Ota fc lye 2012; Shibuya et al. 20lJ). All these Lya sur 



veys were made with ground-based instruments owing to 
their large field-of- views (FOVs) and relatively low sky- 
background in OH-dark windows. 

A complementary way to find high-redshift galax- 
ies is the d ropou t technique (Ste idel fc Hamilton! 1993: 
iGiavaliscol [2002). It has produced a large num- 
ber of Lyman-break galaxies (LBGs) or candidates 
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at z > 6. While large- area ground-based ob- 
servations are efficient to select bright LBGs (e.g., 
Nagao et al.l 120071: iBowler et al.ll2Q12l: ICurtis-Lake et all 



20121 : iHsieh^TaTl 1^)^ the ma- 

jority of the faint LBGs at z > 6 known so far 
were discovered by HST (e.g.. IBunker et all l2QQ4t 
iDickinson et!H 120041 lYan fc Windho rst 2004). In par- 
ticular, with the power of the new HST WFC3 in- 
frared camera, a few hundred gala xies at 6 < z < 12 
have been detected recently (e.g., IBunker et al.l 120101: 
Finkelstein et all 2010: McLureet al. 2010; Oesch et all 
2010; Wilkins et al. 2010; Yan et al. 2010; Bouwens et al. 
201ll : lLorenzoni et aLll2011b lEllis et al.ll2013l ). Although 
most of the LBG candidates are too faint for follow- 
up identification, some of them have been spectro- 
scopically c onfirmed with deep gro u nd-based observa- 
tions (e.g., IStark et al] |2Qlfl UQUi jjiang et all 120111: 



only a few million years (e.g., Pirzka fet al.l 12007). As 
the number of galaxies known at z > 6 increases 
steadily, extensive studies are being carried out on vari- 



Vanzella et al. 2011; Pentericci et al, 



2012; Schenker et al 2012). 



2011} lOno et all 



With the large samples of galaxies (and candi- 
dates) at z > 6, the UV continuum luminosity func- 
tion (LF) and Lya LF are b eing established (e.g. , 
Hu et alJEOlOt iJiarig etlI3|2flll; [K ashikaw a^aD2011l: 
Bradley et all 1201 2: Bouwens et al. 2012a; iHenrv et all 
20121 : lLaporte et al . 1 120121 : IQesch et al.ll2012l ). It is found 



that the faint-end slope of the UV LFs are very steep 
(a ~ —2), so dwarf galaxies could provide enough UV 
photons for reionization (also depending on the ionizing 
photon escape fraction and the ionized gas clumping fac- 
tor). For individual galaxies, their physical properties 
are also being investigated. At z > 6, the rest-frame 
UV/optical light moves to the infrared range. Therefore, 
infrared observations, especially the combination of HST 
and Spitzer Space Telescope (Spitzer), are essential to 
measure properties of these high-redshift galaxies. HST 
near-IR data constrain the slope of the rest-frame UV 
spectrum and help to decipher the properties of young 
stellar populations. Recent HST WFC3 near-IR obser- 
vations have shown that high-redshift LBGs have blue 
rest-frame UV slopes (3 (fx oc \@). The typical val- 
ues of are close to or smaller (bluer) than -2 (e.g., 
Wilk ins et a,l-ll2011b iBouwens et al.| [2Q12bt iDunlop et all 



2012ai iFinkelstei n et al.l I2012D . Though the selection 
criteria of LBGs and the approaches to UV slopes are 
slightly different in the different studies in the litera- 
ture, it is believed that z > 6 galaxies are bluer than 
low-redshift galaxies, implying that they are generally 
younger and have low dust extinction. 

More detailed physical properties (e.g. age and stel- 
lar mass) of high-redshift galaxies have come from SED 
modeling based on broad-band photometry of HST and 
Spitzer as well as optical data. Spitzer IRAC pro- 
vides mid-IR photometry. When combined with HST 
near-IR data, it measures the amplitude of the Balmer 
break and constrains the properties of mature popu- 
lations. The early results on z ~ 6 LBGs snowed 
IRAC detections, and emphasized galaxies wi t h es- 
tablished stellar populations (e.g ., lEgami et al.l |2005; 
lEvles et al.l 120051 : lYan et all 12005ft . Later studies found 
that most LBGs were actually not detected in moder- 
ately deep IRAC images, suggesting that these galax- 
ies were conside rably young er and less massive (e.g., 
lYan et all 12001 lEvles et al.ll2007h . LAEs may repre- 
sent an even more extreme population with ages of 



ous g a laxy samples (e.g., [G onzalez et al. 2010; Ono et aT] 
120101: ISchaerer fc de Barrosl 120101 : iMcLure et all 120111 : 
iPirzkal et al.l [2012b . A diversity of physical properties 
are found among these galaxies, although the Universe is 
younger than one billion years. 

Despite the progress that has been made on proper- 
ties of high-redshift galaxies, there is a lack of a large, 
reliable sample of z > 6 galaxies with spectroscopic red- 
shifts and high-quality infrared photometry. The reason 
is twofold. For galaxies found by ground-based telescopes 
- although they are relatively bright and have spec- 
troscopic redshifts — they often do not have measured 
infrared (rest-frame UV/optical) information. For ex- 
ample, only a small amount of LAEs at z > 6 have been 
observed in the near- IR bands by HST (e.g. JCowie eFap 
l2011aft . On the other hand, galaxies found by HST (e.g. 
those in the GOODS fields), do have near-IR photome- 
try, but most of them are too faint for follow-up spec- 
troscopic identification by current facilities. So they do 
not have secure redshifts, making it difficult for accu- 
rate SED modeling. A model SED derived from a few 
photometric points without a secure redshift is usually 
highly degenerate. Therefore, determination of the phys- 
ical properties of z > 6 galaxies demands a large sample 
of galaxies with both spectroscopic redshifts and high- 
quality infrared photometry. 

In this paper we will present deep HST and Spitzer 
observations of 67 spectroscopically-confirmed (hereafter 
spec-confirmed) galaxies at z > 6. This galaxy sam- 
ple contains 51 LAEs and 16 LBGs, and is the largest 
collection of spec-confirmed galaxies in this redshift 
range. The spectroscopic redshifts of the sample pro- 
vide great advantages for measuring the physical proper- 
ties of galaxies. First, this sample is uncont animated by 
interlopers. A photometrically-selected (hereafter photo- 
selected) sample of high-redshift galaxies may be contam- 
inated by interlopers such as low-redshift red or dusty 
galaxies and Galactic late- type dwarf stars, which may 
cause various bias on measurements of physical proper- 
ties. Second, spectroscopic redshifts remove one criti- 
cal free parameter for SED modeling. For the brightest 
galaxies at z > 6, there are usually only 4-5 broad-band 
photometric points available (e.g., two HST bands, two 
Spitzer IRAC bands, and one optical band) for produc- 
ing stellar population models. Given the limited degrees 
of freedom, a spectroscopic redshift will significantly im- 
prove SED modeling, especially when nebular emission 
is considered. At z > 6, strong nebular emission lines 
such as [Om] A5007, H/3, and Ha enter IRAC 1 and 
2 channels, which affect our estimate of stellar popu- 
lations (e.g., Schaere r fc de B arro si 12010b iFinlator et al.l 



Ifflll : Ide Barros et all 120121 : IStark et all 120131) . Photo- 
metric redshifts with large uncertainties may place these 
nebular lines in wrong bands during SED fitting. With 
spectroscopic redshifts, the positions of nebular lines are 
secured (see Stark et al. 2012 for more discussion). Fi- 
nally, we currently have little knowledge of stellar pop- 
ulations in z > 5.5 LAEs. This is because almost all 
the known LAEs were discovered by ground-based tele- 
scopes, and most of them do not have infrared observa- 
tions. Our paper includes a large sample of LAEs. This 
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allows us, for the first time, to systematically study stel- 
lar populations in LAEs. 

This paper is the first in a series presenting the phys- 
ical properties of these galaxies based on observations 
from HST, Spitzer, and the Subaru telescope. The vast 
majority of the g alaxies were discovered in the Subaru 
Deep Field (SDF; iKashikawa etall l200l . The SDF is 
a unique and ideal field to study the distant Universe: 
it covers an area of over 800 arcmin 2 and has the deep- 
est optical images among all ground-based imaging data. 
The SDF project has been highly successful in search- 
ing for high-redshift galaxies. It h as identified more than 
100 LAEs and LBGs at z > 6 (e.g., Taniguc hTeHIllgO^ : 
Ive et al. 2006; Shimasaku et al. 2006; Nagao et al. 2007; 
Ota et al . 2008; Jiang et al. 2011; Kashika waet al.ll20lH 
Toshikawa et al.ll2012f ). These galaxies are spectroscopi- 
cally confirmed and represent the most luminous galax- 
ies in this redshift range. In this paper we present the 
HST and Spitzer observations, and use these data to- 
gether with the optical data from Subaru to derive the 
rest-frame UV properties, such as their UV continuum 
slopes /3 and Lya luminosities. In the subsequent papers 
of this series we will present rest-frame UV morphology 
and measure stellar populations for these galaxies. 

The structure of the paper is as follows. In Section 2, 
we present our galaxy sample. We also describe our opti- 
cal/infrared observations and data reduction. In Section 
3, we derive the rest-frame UV continuum properties of 
the galaxies. In Section 4, we measure the properties 
of Lya emission lines. We discuss our results in Section 
5 and summarize the paper in Section 6. Throughout 
the paper we use a A-dominated flat cosmology with 
Hp = 70 km s- 1 M pc" 1 , Q m 0.3, and Q A = 0.7 
(Komatsu et al. 2011 ). Al l magnitudes are on the AB 
system (jOke fc Gunnlll983l ). 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. Galaxy Sample 

Table 1 shows the list of the galaxies presented in this 
paper. There are a total of 67 spec-confirmed galaxies at 
5.6 < z < 7: 62 of them are from the SDF and the re- 
maining 5 a re from the Subaru XM M-Newton Deep Sur- 
vey fSXDS: iFurusawa et al1l2008l ) field. They represent 
the most luminous galaxies in terms of Lya luminosity 
for LAEs or UV continuum luminosity for LBGs in this 
redshift range. 

The SDF covers an area of ~ 876 arcmin 2 , and its op- 
tical imaging data were taken with Sub aru Suprime-Cam 
in a s eries of broad and narrow bands (jKashikawa et all 
2004). Especially noteworthy are the deep observations 
with three narrow-band (NB) filters, NB816, NB921, 
and NB973, corresponding to the detection of LAEs at 
z ~ 5.7, 6.5, and 7. The full widths at half maxi- 
mum (FWHM) of the three filters are 120, 132, and 
200 A, respectively. So far the SDF project has spec- 
troscopically confirmed more than 100 LAEs at z c± 
5.7, 6.5, 7, and more than 40 LBGs at 6 < z < 7. 
Our SDF galaxy sample contains 22 LAEs a t z ~ 5.7 
([Shimasaku et alJl2QQ6t IKa shikaw a et al.ll2011ft. 25 LAEs 
at z ~ 6.5 ( Taniguchi et al .1 1200 5: Kashika wa et al.l l2006. 
[20lTh . and the z = 6.96 LAE (jive et al.ll2006l ). The LAEs 
at z ~ 5.7 and 6.5 were selected in a similar way, and 
have a relatively uniform magnitude limit of 26 mag in 



the narrow bands NB816 and NB921, so they are from 
a well-defined flux-limited sample. O ur SDF sample also 
conta i ns 14 LBGs at 5.9 < z < 6.5 ([N agao et aLl l2004 
I2005L 120071 IQta et al.l l2008t Lliang et all 120111 ): two of 
them have not been previously published. These LBGs 
were selected with different criteria and thus have in- 
homogeneous depth. The z'-band magnitud e limit is 
26.1 m ag in Nagao et a l. papers, 26.5 ma g in IQta et al.l 
(|2008[ ). and 27 mag in Uiang et all ([201 lh . For the de- 
tails of candidate selection and follow-up spectroscopy, 
see the original galaxy discovery papers above. 

Five galaxies in our sample , including 2 LBGs at 
z^6jCurtis-La ke et all 1201 2[ ) and 3 LAEs at z ~ 6.5 
(jOuchi et al.ll201Q[ ). are from the SXDS. The two LBGs 
ar e very bright (z f < 25 m ag) and have been analyzed 
bv lCurtis-Lake et al.l (|2013l ). The SXDS optical imaging 
data were also taken with Subaru Suprime-Cam in the 
same broad and narrow b ands, but cover five tim es larger 
area than the SDF does ([Furusawa et al.ll2008[ ). 

In Table 1, the first 62 galaxies are from the SDF and 
the last 5 galaxies are from the SXDS. Columns 2 and 3 
list the object coordinates. They are re-calculated from 
our own stacked optical images (see Section 2.2), and 
are fully consistent with those given in the galaxy dis- 
covery papers (typical difference < 0'/l). Column 4 lists 
the redshifts, measured from the Lya emission lines of 
the galaxies. The galaxy discovery papers used differ- 
ent values for the rest-frame Lya wavelengths (i.e. 1215, 
1215.67, or 1216 A) for redshift measurements. We con- 
vert all these redshifts to the redshifts listed in Table 1 by 
assuming a rest-frame Lya line center of 1215.67 A. The 
last column shows the reference papers of the objects. 
We will explain the rest of the table in the following sub- 
sections. 

2.2. Optical Imaging Data 

The galaxies in this paper were discovered in the SDF 
and the SXDS. The SDF public imaging data include 
stacked images in five broad bands BVRi'z' and two nar- 
row bands NB816 and NB921. The flux limits for these 
data are given in Kashika wa et al.l (|2004f ). However, the 
public data do not include the data taken recently (also 
by Subaru Suprime-Cam), and the data taken in the 
y and NB973 bands. Therefore, we produced our own 
stacked images by including all available images in the 
Subaru archive, as explained below. 

Our data processing began with the raw images ob- 
taine d from the archival server SMOKA (jBaba et al.l 
2002). Raw images with point spread function (PSF) 
sizes greater than l'/2 were rejected. The data were 
reduced, re-sampled, and co-added using a combina- 
tion of t he Suprime-Cam Deep Field REDu ction package 
fSDFRED: lYagi et al.ll200llOuchi et al.ll2004D and our own 
IDL routines. Each image was bias (overscan) corrected 
and flat-fielded using SDFRED. Then bad pixel masks were 
created from flat-field images. Cosmic ray s were identi- 
fied a nd interpolated using L.A. Cosmic (va n Dokkuml 
120011 ). Saturated pixels and bleeding trails were also 
identified and interpolated. For each image, a weight 
mask was generated to include all above mentioned de- 
fective pixels. We then used SDFRED to correct the im- 
age distortion, subtract the sky-background, and mask 
out the pixels affected by the Auto-Guider probe. After 
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individual images were processed, we extra cted sources 
using SExtractor (jBertin fc Arnoutsl Il996[ ) and calcu- 
lated astrome tric and photometric solutions using SCAMP 
(Bertin 2006) for the final image co-addition. In order 
to match the public SDF images, we used the public R- 
band image as the astrometric reference catalog in SCAMP. 
Both science and weight-map images were scaled and up- 
dated using the astrometric and photometric solutions 
measured from SCAMP. Instead of homogenizing PSFs to 
a certain value as was done for the public data, we incor- 
porated PSF information into the weight image so that 
the weight was proportional to the inverse of the square 
of the PSF F WHM. We re- samp led and co-added images 
using SWARP (Ber tin et al.ll2002[ ). The re-sampling inter- 
polations for science and weight images were LANCZOS3 
and BILINEAR, respectively. We modified SWARP so that 
it performs sigma-clipping (5a rejection) of outliers when 
co- adding images. 

We produced images in the six broad bands BVRi'z'y 
and three narrow bands NB816, NB921, and NB973. The 
final products are a stacked science image, and its cor- 
responding weight image, for each band. The pixel and 
image sizes of our products are the same as those of the 
public SDF images. The typical PSF FWHM is / /6-0 / /7. 
We determined photometric zero points from the public 
SDF images whose zero points were known. The SXDS 
optical imaging data were reduced in the same way as 
we did for the SDF. 

Our SDF broad-band images are deeper than the pub- 
lic images in the Ri V bands, because the SDF team 
obtained more imaging data from Subaru recently, and 
we included all the available data above. For example, 
the public z'-band image has a depth of z' = 26.1 mag 
(5a detection in a 2" diameter aperture) with a total in- 
tegration time of ~7 hr. To date, the total z'-band data 
amount to ~29 hr. Our new z'-band image has a depth 
of ~27.0 mag, 0.9 mag deeper than the public image. In 
addition, the public data do not have images in the y 
and NB973 bands. The SDF ?/-band imaging data were 
taken with two different sets of CCDs (MIT and Hama- 
matsu), and have been used to search for z ~ 7 LBGs 
(jOuchi et all 12009( 1. Although the fully-depleted Hama- 
matsu CCDs are twice as sensitive as the MIT CCDs in 
the y band, their sensitivity (quantum efficiency) curves 
have similar shapes. Therefore, w e stacked all t he im ages 
taken with both sets of CCDs, as lOuchi et al.l (|2009f ) did 
in their work. The total integration time was 24 hr, and 
the depth was about 26.0 mag. 

In Table 1, Columns 5 to 7 show the total magni- 
tudes of the galaxies in one of the three narrow bands 
(NB816, NB921, or NB973, depending on redshift) and 
two broad bands z'y. For a given object in a given fil- 
ter, its total magnitude is measured as follows. We first 
use SExtractor to measure the aperture photometry of 
the object in a circular 2 / /0 diameter, after the local sky- 
background is subtracted. In cases where galaxies are 
close to their neighbors, we use a circular l'/6 diameter. 
Then an aperture correction is applied to correct for light 
loss. The aperture correction is determined from more 
than 100 bright, but unsaturated point sources in the 
same image. Almost all our galaxies are point-like ob- 
jects and are not resolved in the SDF images. In cases 
where galaxies clearly show extended features, we use 
MAG_AUTO magnitudes as the best estimates of the 



total magnitudes. In Table 1, we did not list broad-band 
photometry for the galaxies with weak detections (< 5<r) 
in the J band. These galaxies are excluded in most of 
our analysis in this paper. Detections below 3a in any 
of these bands are not listed as well. We carefully check 
each object, and find that object no. 12, a z ~ 5.7 LAE, 
is overlapped by a bright foreground star (this was also 
pointed out in the galaxy discovery paper). The separa- 
tion between the two objects is only ~ 0'/25, so we are 
not able to do photometry for this galaxy. We will not 
use this galaxy in the following analysis. 

Figure 1 shows the transmission curves of six Subaru 
Suprime-Cam filters, including three broad-band filters 
i' z'y and three narrow-band filters NB816, NB921, and 
NB973. We do not show the other three bluer filters 
BVR, because our galaxies were not detected in these 
bands. In Appendix A, we show the thumbnail images 
of the galaxies in one of the three Subaru Suprime-Cam 
narrow bands and the three broad bands z, 2/, and y. 

2.3. HST Near-IR Imaging Data 

We obtained HST near-IR imaging data for a large 
sample of spectroscopically confirmed SDF galaxies at 
5.6 < z < 7 from three HST GO programs 11149 (PI: 
E. Egami), 12329 and 12616 (PI: L. Jiang). The three 
programs together with our Spitzer programs (see the 
next subsection) were designed to characterize the phys- 
ical properties of these high-redshift galaxies. The HST 
observations were made with a mix of strategies. As 
the proposal for GO program 11149 was written, HST 
WFC3 had not yet been installed. Our plan was to ob- 
serve 20 bright galaxies with NICMOS in the F110W 
(hereafter Jno) and F160W (hereafter #160) bands. The 
observations were made with camera NIC3. The expo- 
sure time was two orbits per band per pointing (or per 
object), split into four exposures. The total on-source 
integration per band was about 5400 sec. We performed 
sub-pixel dithering during the four NICMOS exposures 
to improve the PSF sampling and to remove bad pixels 
and cosmic rays. 

After we finished several targets, NICMOS failed, and 
then WFC3 was installed in May 2009. Since WFC3 
has much higher sensitivity, the rest of the galaxies in 
GO 11149 and all the galaxies in GO 12329 and 12616 
were observed with WFC3. Because of the larger FOV 
of WFC3 and the increasing number of known z > 6 
galaxies in the SDF, we changed our observing strategy 
for GO 12329 and 12616. We chose the regions that cover 
at least 3 galaxies (up to 5) per HST /WFC3 pointing. 
By doing this we were able to observe many galaxies in 
a small number of telescope pointings. We used F125W 
(hereafter J125) and Hiqq in these two programs. In GO 
12329, the exposure time was one orbit per filter per 
pointing. We found that some galaxies in our sample 
were barely detected with this depth. Therefore in GO 
12616, all new WFC3 pointings had a depth of two orbits 
(on-source integration ~ 4800 sec; see Windhorst et al. 
2011 for the anticipated sensitivity in two-orbit WFC3 
exposures). We also re-observed most galaxies that were 
observed in GO 12329, so that they have a total of two- 
orbit depth per band. As in the NICMOS observations 
above, the WFC3 observation of each pointing (per band) 
was split into four exposures. Sub-pixel dithering during 
the exposures was also performed to improve the PSF 
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sampling and to remove bad pixels and cosmic rays. The 
transmission curves of the three WFC3 filters are plotted 
in Figure 1. 

As we mentioned above, several galaxies in our sam- 
ple were found in the SXDS. They were covered by the 
UKIDSS Ultra-Deep Survey (UDS). Their HST WFC3 
near-IR data were obtained from the Cosmic Assembly 
Ne ar- i nfrared Deep Extragalactic Legacy Survey ( C AN- 
DELS: iGrogin et alJl2QlH : iKoekemoer et al.ll2011h . The 
exposure depth of the CANDELS UDS data is 1900 sec 
in the J125 band and 3300 sec in the Hiqq band. They 
are shallower than our WFC3 data for the SDF. 

In Appendix A we show the thumbnail images of the 
galaxies in the two HST near-IR bands Jno (or J125) and 
Hiqq. Note that the Ji25-band image of the z = 6.96 
LAE (no. 62) was obtained from GO program 11587 
(jCai et al.ll20lTl ). This image also covers another two 
objects (no. 2 and no. 59) in our sample. Due to the 
high sensitivity of WFC3, the majority of the galaxies 
in our sample are detected with high significance in the 
near-IR images. Only 15 of them — among the faintest 
in the optical bands — have weak detections (< 5a) in 
the J band. 

Our WFC3 data processing began with individual flat- 
fielded, flux calibrated WFC3 infrared exposures deliv- 
ered by the HST archive. Each pointing used a standard 
4-point dither pattern to obtain uniform half-pixel sam- 
pling. All exposures for each poi nting were combined us - 
ing the software Multidrizzle (Koek emoer et al. 2002), 
with an output plate scale of / /06 per pixel and a 
pixfrac parameter of 0.8. For our observations, this 
provides Nyquist sampling of the PSF and relatively uni- 
form weighting of individual pixels. RMS maps were also 
derived from the Multidriz zle weight- map i mage s, fol- 
lowing the procedure used in lDickinson et al] ((2004) that 
accounts for correlated noise in the final images. 

Our NICMOS NIC3 data w ere reduced using a fully 
automated pipeline NICRED (jMagee et al.ll2007D . NI- 
CRED starts its processing with raw science data and 
calibration files, and sequentially runs a set of steps 
from basic calibration to final mosaic generation. The 
most important part of the reduction is the basic cali- 
bration, because the NICMOS data suffer from variable 
bias, electronic ghosts, and cosmic-ray persistence. The 
basic procedure is as follows. NICRED first removes elec- 
tronic ghosts, handles variable bias, rejects cosmic rays, 
and subtracts sky background. It then removes possi- 
ble cosmic-ray persistence and residual instrument sig- 
natures. The rest of the procedure are common steps, 
including flat fielding, a non-linearity correction, image 
alignment and registration, etc. It finally creates driz- 
zled images and the corresponding weight maps. The 
pixel size in the final products is 0'/l, a half of the native 
size. 

We list the near-IR photometry of the galaxies in 
Columns 8 to 10 in Table 1. All galaxies were observed 
in the Hiqo band and one or two of the J bands (Jno 
and/or J125). We perform photometry using SExtractor 
in dual-image mode with RMS map weighting. The J- 
band is used as the detection image. A single set of 
parameters is used for detection in all images. We com- 
pute the total magnitudes of the objects in both filters 
by using MAG_AUTO elliptical apertures, with a Kron 
factor of 1.8 and a minimum aperture of 2.5 semi-major 



radius, then applying aperture corrections. The aperture 
corrections are calculated by running SExtractor again 
with five Kron factors between 2.0 and 4.0, analyzing 
thousands of sources in the range J < 26 mag. We cal- 
culate a sequence of median correction vs. Kron factor, 
which we extend to infinite radius with an exponential 
fit. All final photometry apertures are visually inspected, 
to screen against contamination from nearby objects or 
spurious detections. Detections below 5<r in the J band, 
or detections below 3a in the H band, are not listed in 
Table 1. 

2.4. Spitzer Mid-IR Imaging Data 

We obtained Spitzer IRAC imaging data for the SDF 
from two GO programs 40026 (PI: E. Egami) and 70094 
(PI: L. Jiang). In GO 40026, we observed 20 luminous 
galaxies at z > 6, which were also the targets of our 
HST program 11149. The observations were made in 
IRAC channel 1 (3.6 /im) and channel 2 (4.5 /im). The 
other two channels (5.8 and 8.0 /im) were observed si- 
multaneously. We used the cycling dither pattern with 
medium steps and a frame exposure time of 100 sec. The 
on-source integration time per channel per pointing is 3 
hr. In GO 70094 (during the Spitzer warm mission), 
we modified our observing strategy and observed high 
galaxy surface density regions, so that we can observe a 
large number of galaxies with slightly more than 10 tele- 
scope pointings. The cycling dither pattern with small 
steps was used with a frame exposure time of 100 sec. 
We focused on IRAC channel 1 and increased on-source 
integration time to ~6 hr per pointing. The two pro- 
grams imaged roughly 70% of the SDF to a depth of 3-6 
hours. They covered all the 62 SDF galaxies in channel 
1 and 51 (out of 62) galaxies in channel 2. 

The IRAC Data processing began with the cBCD (cor- 
rected BCD) products and associated mask and uncer- 
tainty images delivered by the Spitzer Science Center 
(SSC) archive. The image fits headers were updated to 
match the astrometry in the SDF optical images. Then 
the images were processed, drizzled, and co-added us- 
ing using the SSC pipeline M0PEX. This is a fully auto- 
mated procedure. The final images have a pixel size of 
0'/6, roughly a half of the IRAC native pixel scale. The 
IRAC images for the SXDS galaxies were obtained from 
the SSC archive (programs 60022 and 80057). The data 
were reduced in the same way above. In Appendix (the 
last two columns) we show the thumbnail images of the 
galaxies in the two IRAC channels. 

Mid-IR photometry is complicated by source confu- 
sion. Galaxies in deep IRAC images are often blended 
with nearby neighbors, so accurate photometry requires 
proper deblending and removing neighbors. We will 
present IRAC photometry of our galaxies in the second 
paper of this series. 

3. REST-FRAME UV CONTINUUM PROPERTIES 

With the narrow- and broad-band photometry from 
the NB816 to the Hiqq band, we are able to measure the 
spectral properties, such as Lya flux and UV continuum 
slopes, for most of the galaxies in this sample. In order 
to measure the rest-frame UV properties, we produce a 
model spectrum for each of these galaxies. The model 
spectrum f ga \ (in units of erg s _1 cm -2 A -1 , like com- 
monly used f\) is the sum of a Lya emission line and a 
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power-law UV continuum with a slope /?, 

f ga i = AxS Lya +Bx \P, (1) 

where A and B are scaling factors in units of 
erg s _1 cm -2 A -1 , and Stye* (dimensionless) is the Lya 
emission line profile. In the redshift range and wave- 
length range considered here, other neb ular lines such a s 
Hell A1640 can be safely ignored (e.g., iCai et alJl20li[ ). 
We first determine the continuum parameters B and /? by 
fitting the continuum /con, a = B x (or J COI1)I/ ex A^ +2 ) 
with the broad-band data that do not cover the Lya 
emission. Because AB magnitude tuab and Log(A) have 
the simplest linear relation tuab oc (/3 + 2) x Log(A), the 
continuum fitting is performed on tuab and Log(A) using 
the IDL routine linfit.pro. This routine takes Log(A), 
mAB, and the errors in mAB, and performs a standard 
(weighted) least-squares linear regression by minimizing 
the chi-square error statistic. The errors propagate to 
the final parameters. 

In order to reliably measure continuum slopes, we only 
consider the galaxies with > 5a detections in the J band 
(the deepest band). We further require that these galax- 
ies have > 3cr detections in at least one more broad band 
that does not cover the Lya emission. Detections be- 
low 3<j are not used here. For the galaxies at z < 6.1, 
four broad bands z'yJH can be used. As shown in Table 
1, the z < 6.1 galaxies with > 5a detections in the J 
band and with > 3a detections in another band are usu- 
ally detected (> 3a) in the other two bands, so all four 
broad bands are used for most of the galaxies at z < 6.1. 
For the galaxies at z > 6.1, the z' band covers Lya, so 
only three broad bands yJH can be used. However, al- 
most half of the z > 6.1 galaxies with J-band detections 
(> 5a) are not detected at > 3a in the y band. For 
these galaxies, we use only the JH bands to derive their 
UV slopes. The information (or upper limits) from the 
?/-band photometry is ignored, because the Ji^-band im- 
ages are much deeper (by 1-1.5 mag) than the y-band 
images, and the ^/-band upper limits provide little addi- 
tional constraint on /?. In addition, in rare cases where 
some galaxies in the crowded optical images are signifi- 
cantly affected by their environments like bright neigh- 
bors (e.g., no. 2 and no. 4), their optical data are not 
used. 

We then determine parameter A for Lya emission line. 
iKashikawa et al.l (|2Qllh generated two composite Lya 
emission line profiles for their z ~ 5.7 and 6.5 LAEs, 
and found that the two profiles were very similar, show- 
ing no significa nt evolution from z c± 5.7 to 6.5 (see also 
iHu et al.ll2010f ) . We use the two composite line profiles as 
our model Lya emission lines in Eq. 1: we use the z ~ 5.7 
profile for the galaxies at z < 6.2 and the z ~ 6.5 profile 
for the galaxies at z > 6.2. The shape of the Lya emis- 
sion lines has negligible impact on our results, because 
even the three narrow-band filters used for Lya detection 
are much wider than the typical Lya line width. We ap- 
ply IGM absorption to the model spectrum (the IGM 
absorption has already been taken into account in the 
composite Lya emission lines). The neutral IGM fraction 
increases dramatically from z = 5.5 to z = 6.5, causing 
complete Gunn-Peterson troughs in some z > 6 quasar 
spectra (jFan et al.ll2006[ ). It is thus important to include 
IGM absorption in the model spectrum. We calculate 



IGM absorption in the same way as iFan et al.l J2001). 
The IGM-absorbed model spectrum is convolved with 
the total system response, such as filter transmission and 
CCD quantum efficiency. Finally, the Lya emission flux 
(or parameter A) is calculated by matching the model 
spectrum to the narrow-band photometry (for LAEs) or 
the z'-band photometry (for LBGs). For LBGs, we will 
see in Section 4.1 that even small z'-band photometric er- 
rors will cause large uncertainties on measurements of A, 
so we adopt the Lya flux of the LBGs from the galaxies 
discovery papers, which were derived from deep optical 
spectra. 

Figure 2 illustrates our procedure by showing our 
model fit to the photometric data points of two LAEs 
at z ~ 5.7 and z c± 6.5, respectively. For the z ~ 5.7 
LAE in this specific case, we use the photometric data 
in the z'yJuoHiQo bands (red circles) to fit the model 
continuum. When the continuum (strength B and slope 
j3) is derived, we apply IGM absorption. The continuum 
flux at the blue side of the Lya line (the dashed line 
in the figure) is absorbed. Finally the factor A is de- 
termined by scaling the Lya emission line to match the 
NB816-band photometry. For the z c± 6.5 LAE in Figure 
2, we use the data in the yJi25Hi60 bands to compute 
its continuum, and use the NB921-band photometry to 
measure its Lya emission. After A, £?, and j3 are de- 
termined, other physical quantities such as Lya and UV 
luminosities, rest-frame Lya equivalent width (EW), and 
star formation rates (SFRs) are also calculated. The re- 
sults are shown in Table 2. 

3.1. UV Continuum Slope /3 

The rest-frame UV-continuum slope provides key infor- 
mation to constrain young stellar populations in galaxies. 
For z > 6 LBGs in the literature, especially those se- 
lected from HST photometric samples, their UV slopes 
are usually measured from two broad bands J125 and 
Hiqo that do not cover Lya emission. Recent stud- 
ies show that high-redshift LBGs general ly have steep 
UV slopes and low dust content (e.g., [B ouwens et al.l 
[2QQ9ll2Q12lJ _ lDunloD e t al. 2012a; Finkelstein etaIll2Ql3 
IGonzalez et al.l 12012): IWalter et al.l l2012f ). For exam- 
ple. iBouwens et al.l (j2009[ ) found that 3 decreases from 
3 ~ — 1-5 at z ~ 2 to /3 ~ — 2 at z ~ 6. Finkels tein et al.l 
(|2Q12f ) showed that /3 changed from /3~-1.8at2~4 
to () ~ —2.4 at z ~ 7. Extremely steep slopes of 
/3 ~ — 3 have been reported i n very faint z ~ 7 galax- 
ies (e.g., Bo uwens et al. 2010). It was also found that 
lower - luminosity galaxies tend to have steeper slopes 
(e.g., IBouwens et al.l 120121 IGonzalez et al.ll2012[ ). Re- 
cently the results on the bl ue colors in high-red shift 
LBGs have be e n ques tioned bv lDunlop et al.l (j2012af ) and 
iMcLure et al.l (|2011l ). who claimed that previous mea- 
surements were significantly affected by contaminants 
in the photo-selected LBG samples and the low signal- 
to- noise (S/N) ratios of faint objects, and that the ex- 
tremely steep slopes were caused by the so-called "/? 
bias". This bias is partly produced by the removal of 
galaxies with larger (redder) /3 values from the sample as 
lower-redshift interlopers, which would preferentially re- 
move flux-boosted sources in the -H160 band and skew the 
/3 distribution to smaller (i.e., bluer) v alues as a result. 
As pointed out bv IBouwens et all (|2012bl ). this /? bias can 
be large if the same bands (J125 and Hiqq) are also used 
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to select LBGs. With a more restricted sele ction of LBGs 
candida tes from the same HST data set, iDunlop et al.l 
(|2012a[ ) and lMcLure et al.l ([201 ID found that the average 
slope of z ~ 7 LBGs is /3 ~ —2, which is not steeper 
than the bluest galaxies at z < 4. They further found 
no relation between /? and the UV contin uum luminos- 
ity. A s the debate is sti ll going on (e.g., IB ouwens et al.l 
I2012bl : iFinkelstein et al.ll2Q12l : iDunlop et aHl2Q12b[ ). the 
key is to construct a high-redshift galaxy sample that is 
free from such a bias. 

Our galaxy sample presented here is unlikely to be 
strongly affected by this /3 bias, not only because these 
galaxies have spectroscopic redshifts but also because our 
original source selection did not utilize any information 
on the rest-frame UV continuum slope. Furthermore, 
we did not include weak detections in our analysis. For 
many galaxies, especially the z ~ 5.7 LAEs and z ~ 6 
LBGs, f3 was measured from more than two bands. It 
should be pointed out, however, our LBGs were spec- 
confirmed to have strong Lya emission, which could bias 
the LBG sample to galaxies with lower extinction and 
thus to those with bluer UV continuum slopes. 

Figure 3(a) shows the rest-frame UV continuum slopes 
f3 as a function of M1500, the absolute AB magnitude of 
the continuum at rest-frame 1500 A. M1500 is calculated 
from B x in Eq. 1. Figure 3 only includes the galaxies 
that have measured /3 in Table 2. These galaxies cover 
the brightest UV luminosity range M1500 < —19.5 mag. 
For the galaxies that do not have measured /?, they do 
not have reliable detections (< 5<r) in the near-IR bands, 
so they are usually fainter than M1500 = —19.5 mag (see 
Section 4.2 and Figure 6). The green circles in Figure 
3(a) represent the LBGs at z ~ 6, and the blue and red 
circles represent the LAEs at z ~ 5.7 and 6.5 (including 
z ~ 7), respectively. The values of Mi 500 and j3 are also 
listed in Columns 2 and 3 of Table 2. Our measurements 
of f3 are reliable because of the secure redshifts and multi- 
band photometry. As we addressed above, the slopes are 
calculated from all available broad-band data that do not 
cover the Lya emission line. For any f3 that was derived 
from more than two bands, its error is usually smaller 
than 0.4. So most of the z ~ 5.7 LAEs and z ~ 6 LBGs 
have relatively small errors on j3. For the z ~ 6.5 LAEs, 
their slope errors are generally larger, as many of them 
were only detected in the two near-IR bands in addition 
to NB921. 

The galaxies in Figure 3 apparently have steep 
UV-continuum slopes that are roughly between /3 ~ 
— 1.5 and /3 ~ —3.5. The weighted mean is /3 = 
— 2.29, and the median value is j3 = —2.35. This 
slope is slightly steeper than those from previous ob- 
servat ion s (e.g., IBouwens et al.l l2012bl : IDunlop et al.l 
2012allbl: I Finkelst ein et al.l 12012ft and simul ations (e.g., 



Finlator et al.l 2011; Dayal & Ferrara 2012) in the lit- 



erature. The slope f3 also appears to show weak de- 
pendence on Mi 500 that lower-luminosity galaxies have 
steeper slopes. The dashed line is the best linear fit to 
all the data points. This trend is caused by the several 
most luminous galaxies at M1500 — —22 mag. These lu- 
minous galaxies have red colors /3 ~ — 1.8. If we exclude 
these galaxies in our fit, the trend almost disappears. 
The dash-dotted line in the figure illustrate the result by 
displaying the best fit to the galaxies at M1500 > —21.7 



mag, and its slope is consistent with zero. 

In the above analysis, galaxies with near-IR detections 
below 5<j in the J band were discarded. This has negligi- 
ble impact on our results. Our HST near-IR data have 
relatively uniform depth (two orbits per band per point- 
ing), so the 5cr detection cut indeed puts a flux limit on 
M1500 in Figure 3(a), which does not affect our results 
in the bright region Mi 500 < —19.5 mag. We also per- 
form two simple tests to investigate whether large errors 
on photometry or f3 have strong impact on our results. 
In the first test we stack the images of 9 galaxies whose 
errors of f3 (a (3) are greater than 0.8 in Figure 3, and 
then measure photometry and calculate j3 on the stacked 
images in the J125 and Hiqo bands. The input images 
are scaled so that the galaxies have the same magnitude 
in the J125 band before stacking, otherwise the final im- 
age would be dominated by the brightest objects. The 
stacked images have much higher signal-to- noise (S/N) 
ratios. The 9 galaxies have a median slope of j3 c± —2.1. 
The slope for the stacked object is /3 = —2.2 ±0.2, which 
agrees with the median value of the individual galaxies. 
The result is shown as the magenta star in Figure 4, 
where the gray circles represent the same galaxies shown 
in Figure 3(a). In the second test there are 10 galaxies 
with j3 < —2.8 that have both J125 and Hiqo images. The 
median value of their slopes is j3 c± —3.0. We stack their 
images in the same way as we did above. The slope for 
the combined object is /3 = — 2. 9 ±0.2, which is also con- 
sistent with the median value of the individual galaxies. 
The magenta square in Figure 4 represents the combined 
object. From these tests, our measurements of j3 are not 
biased even for faint objects with relatively large photo- 
metric errors. 

3.2. /3-M1500 Relation 

It is generally believed that high-redshift galaxies have 
steep UV slopes /?, and j3 changes with redshift: higher- 
redshift galaxies tend to be bluer due to less dust ex- 
tinction. The relation between j3 and UV luminos- 
ity (Muv, or M1500 in this paper) is still controver- 
sial. In Fi gure 4 we compare our results with re- 
cent studies (IBouwens et al.ll2012b1: IDunlop et al.ll2Q12at 
IFinkelstein et al.ll2012f ). IBouwens et al.l (j2012bf )found a 
correlation between (3 and Muv at all redshifts in a range 
of 4 < z < 7. They claim that galaxies have bluer UV 
colors at lower luminosities. The green triangles in Fig- 
ure 4 represent the average slopes at z ~ 6 from their 
study. IDunlop et al.l (|2012al ) used a more restricted se- 
lection of LBGs by excluding low S/N detections from 
a similar HST data set. They found that high-redshift 
galaxies have an average slope of j3 ~ —2, and it does 
not change with Muv- The blue crosses in Figure 4 rep- 
resent the averag e slope s at 5 < z < 7 from their study. 
IFinkelstein et all (|2012f ) also found that j3 shows minor 
relation with M1500 in the redshift range from z = 4 to 
8. The red diamonds in Figure 4 represent the average 
slopes at z c± 6 from their study. 

Figures 3 and 4 suggest that in our sample there is little 
dependence of f3 on M1500 in the range of M1500 < —19.5 
mag. The three studies mentioned above were all based 
on HST data that are ~ 1 — 2 mag deeper than ours. 
In the luminosity range of M1500 < —19.5 mag, these 
studies do not show evidence of a significant correlation 
either, though there is likely a weak trend at the fainter 
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range M1500 > —19.5 mag. Therefore, it is possible that 
the reported correlation between j3 and Muv, if exists, 
is mainly driven by very faint galaxies. 

The mean value of our measur ed /3 is -2.3. It is slightly 
steeper than /3 = 2.1 - 2.2 by iBouwens et~aT1 (|2Q12bf ) 
in the sam e luminosity range, an d is also steeper than 
~ -2 by iDunlop et all (|2012af ) and iFinkelstein et al.l 
(2012). This is due to the nature of our sample. The 
galaxies in the previous studies are photo-selected LBGs. 
Our galaxies are spectroscopically confirmed and have 
strong Lya emission lines, suggesting that they likely 
have lower dust extinction, and thus steeper UV slopes. 

3.3. Slopes in LAEs 

It is not clear whether high-redshift LAEs and LBGs 
represent physically different galaxy populations. Early 
observations found that LAEs at z = 5 ~ 6 are smaller 
and less massive compared to LBG s, and thus constitute 
younger stellar populati ons (e.g., iDow-Hvgelund et al.l 
120071 : iPirzkal et aH 120071 ). However, this apparent dif- 
ference could be a direct result of selection effects, and 
LAEs may represen t a subset of LBGs wit h strong Lya 
emission lines (e.g. JDaval fc Ferra ra 2012). Whether or 
not LAEs and LBGs have similar populations will be 
reflected in their physical properties. For example, if 
LAEs are substantially younger, they are expected to 
have steeper UV slopes. 

In Figures 3 and 4 the distribution of our slopes for 
the z c± 5.7 and z ~ 6.5 LAEs are = —2.3 ± 0.5 and 
= — 2.0 ± 0.8, respectively. The distribution for the 
LBGs is /3 = —2.5 ±0.4. Although the scatters are large, 
the LAE slopes are apparently not steeper than the LBG 
slopes in our sample. They are also not steeper than 
the slopes of the LBGs in previous studies mentioned 
above, which is contrary to the assumption that LAEs 
contain younger stellar populations than LBGs do. This 
indicates that most LAEs are probably not as young as 
what was previously suggested. 

3.4. Galaxies with Extreme UV Slopes 

The majority of our galaxies have slopes > —2.7 in 
Figure 3(a). There are a fraction of galaxies whose slopes 
are steeper than -2.7. Some of them have relatively small 
errors around 073 ~ 0.4. The existence of these ~ — 3 
galaxies are statistically robust (not a distribution tail). 
In Figure 3(b) our statistical test shows that the distribu- 
tion of the observed is broader and flatter than the dis- 
tribution of expected if all the galaxies have an intrinsic 
= —2.3 with the observed uncertainties. It clearly in- 
dicates a statistically significant excess of galaxies with 
/3 ~ — 3. These galaxies do not show distinct properties 
in many aspects such as Lya EW and UV morphology. 
The existence of these c± — 3 galaxies in our sample 
is intriguing. The UV colors cannot be arbitrarily blue 
for a given initial mass function (IMF). For the com- 
monly used Salpeter IMF, a /3 ~ — 3 usually means a 
very young stellar population with extremely low metal- 
licity and no dust content, regard l ess of star- formatio n 
histories (e.g., iFinlator et aLl[2QlTt iWilkins et al.ll2011[ ). 
In particular, the < — 3 galaxies, if confirmed by fu- 
ture deeper observations, will be challenging for current 
simulations (jFinlator et al.ll2"01l[ ). In addition, the broad 
distribution of the observed also suggests significant in- 
trinsic scatter in dust and/or age. Detailed discussion is 



deferred to the next paper when we perform SED mod- 
eling with the IRAC data. 

From Figure 3 the five most luminous galaxies at 
M1500 — —22 mag have UV slopes ~ —1.8, which is 
significantly redder than the mean value = —2.3 of the 
sample. The measurements of are robust since these 
galaxies are very bright and detected in almost all the 
broad bands redward of Lya. In particular, they are all 
detected at high significance in the IRAC bands. Their 
Lya emission is relatively small compared to their UV 
continuum emission. From the HST images, they have 
very extended morphologies, even with double or mul- 
tiple cores, suggesting that they are interacting systems 
with enough amounts of dust to redden their UV col- 
ors. These galaxies are similar to t he lum inous galaxies 
at z ~ 6 foun d bv IWillott et all rf2013h in the CFHT 
Legacy Survey. IWillott et al.l (|2013f ) selected a sample of 
luminous galaxies with Muv — —22 in a 4 sq. degree 
area and found that their galaxies have very red colors 
with a typical slope = —1.4, which is caused by large 
dust reddening of Ay > 0.5. 

4. Lya PROPERTIES 

4.1. Lya Luminosity 

Spectroscopic data generally provide reliable line lumi- 
nosities, but this becomes complicated for z > 6 galaxies. 
During the spectroscopic observations of these galaxies, 
the total integration time is at least a few hours per ob- 
ject or per slit mask. Sometimes the images of one object 
are taken in different nights. Light loss varies with many 
parameters, such as varying seeing, different observing 
conditions, possible offsets between targets and slits, and 
even the sizes of targets. So it is difficult to accurately 
correct for light loss. In addition, even with a few hour 
integration on the largest telescopes, the spectral quality 
of many galaxies are still poor, simply because they are 
faint. All of the above effects can easily result in large 
uncertainties in the measurements of their Lya luminosi- 
ties. 

Photometric data provide an alternative way to mea- 
sure the observed Lya flux, especially for LAEs. We 
have high-quality narrow-band and broad-band photom- 
etry and secure redshifts. Our procedure used all avail- 
able data and was straightforward. The Lya luminosities 
of the LAEs were derived from narrow-band photometry 
after their underlying continua were subtracted. For the 
majority of the LAEs in our sample (EW > 20 A), their 
narrow-band photometry is dominated by Lya emission, 
meaning that the real Lya emission is smaller — but not 
much smaller — than the total narrow-band photometry. 
Therefore, we neither significantly overestimate the Lya 
luminosity (the maximum is the total narrow-band emis- 
sion), nor significantly underestimate it (the minimum is 
at least a half of the maximum for most LAEs). 

The measurements of our Lya luminosities for LBGs 
are less accurate. For most of the LAEs, their narrow- 
band photometry is dominated by the Lya emission, so 
their Lya luminosities can be well determined. For the 
LBGs, Lya luminosities were calculated from the z'-band 
photometry. The z' band is very broad, and its photom- 
etry is almost completely dominated by the UV contin- 
uum flux. Therefore, even small uncertainties on the z'- 
band photometry or on measurements of UV continuum 
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will cause very large uncertainties on the Lya emission. 
Therefore for the Lya luminosities of our LBGs, we di- 
rectly adopt the values from the galaxy discovery papers, 
which were derived from optical spectra. All the mea- 
sured Lya luminosities are given in Column 4 of Table 
2. 

Figure 5 compare s our Lya lumi n osities of the LAEs 
to those given in iKashikawa et al.l ([201 lh . The upper 
panel shows the comparison between our luminosities 
and the luminosities from the photometric data in the 
galaxy discovery papers. The photometric data used in 
these papers are usually the narrow-band and z'-band 
photometry. Our results are well consistent with these 
previous measurements with a small scatter. The lower 
panel shows the comparison between our measurements 
and the measurements from the spectroscopic data in 
IKashikawa et al.l (|2011f ). The two results also agree with 
each other, although our measurements are slightly larger 
and the overall scatter is larger than that in the upper 
panel. 

4.2. Lya EW 

Column 5 of Table 2 lists the Lya EWs of the galaxies 
that have measured UV continuum properties in Section 
3.1. These galaxies are also shown as filled circles in Fig- 
ure 6. The bluest galaxy with an unusual f3 = —4.37 is 
excluded. We also estimate EWs for the LAEs that have 
weak or none detections (< 5a) in the J band. These 
galaxies could potentially have very large EWs. In or- 
der to estimate their EWs (or the lower limits of EWs), 
we assume a UV slope f3 = —2.3 and compute UV flux 
from their J-band photometry. Almost all these galaxies 
are detected at > 2a in the J band. If a detection is 
below 2<j, we use 2a as an approximation. The measure- 
ments are shown as the open circles in Figure 6. The up- 
per panel shows the relation between the Lya luminosity 
and UV continuum luminosity M1500. The Lya luminos- 
ity has a weak dependence on M1500. The lower panel 
displays the measured Lya EWs as a function of M1500. 
Our galaxies have moderately strong EWs in the range 
from -10 to -200 A. The median E W value of LAEs 
i s 81 A, consistent with those given in Kashik awa et al.l 
(|2011f ). As expected, the LAEs with weak near-IR detec- 
tions have much higher EWs roughly between 100 and 
300 A. 

Figure 6 (the lower panel) shows that lower-luminosity 
galaxies tend to have higher EWs. It has been a long- 
standing debate whether this is a real relation, or just 
the result of a selection effect, i.e., in a flux-limited sur- 
vey, lower-luminosity galaxies naturally have larger EWs 
in order to be detected in imaging data or identified in 
spectroscopic data. Many papers have reported this in- 
verse relation between Lya EW and UV luminosity to 



be real at z > 3 (e.g., IB haplev et al.l 
2006; Reddy et al. 2008; Stark et al. 
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lAndo et al.l 
while oth- 
ers claimed that the relation is not obvious, when the 
selection effects ar e properly taken into account (e.g., 
iNilsson et all 120091 : ICiardullo et~aT1 120 121 ). In Figure 6 
there is an apparent lack of large-EW LAEs at the high- 
lumin osity end that do es not suffer from this selection 
effect ([Ando et al.l 120061 ). suggesting that the relation in 
our sample is likely real. However, the relation should 
be weaker than it appears, because we have missed faint 



galaxies with small Lya EWs by selection. We will fur- 
ther assess the effects of the source selection on the shape 
of the EW(Lya)-Muv relation in Section 5.1. As we will 
mention in the next subsection, one explanation for the 
relation between Lya EW and UV luminosity (if exists) is 
that the escape fractions of Lya photons in higher-SFR 
(i.e., higher UV luminosity) galaxies are sm aller (e.g., 
iForero-Romero et al.ll20ll IGarel et al.ll2012l ). 

4.3. UV and Lya SFRs 

We estimate SFRs from the UV continuum luminosity 
and Lya luminosity, respectively. The UV SFR is derived 
by: 

SFR(UV) = 1.4 x 1(T 28 L„(UV) M Q yr" 1 , (2) 

where L^fUV) is the UV cont inuum luminosity 
(|Kennicuttl H998: Ma dau et all 11998ft . We use the av- 
erage of the luminosities at rest-frame 1500 and 3000 A 
to compute L zy (UV). For the Lya SFR, we use the ex- 
pression: 

(3) 



SFR(Lya) = 9.1 x 10~ 4d L(Lya) M yr~ J 

which is based on the li ne emission ratio of Lya to Ha in 
Case B recombination ([Osterbrock et all 119891) and the 
relati on between SFR and the Ha luminosity ([Kennicuttl 
1998). The derived SFRs are listed in Columns 6 and 7 
of Table 2. 

Figure 7 compares the two estimated SFRs without 
correction for dust extinction. As in Figure 6, the 
filled circles are the galaxies with reliable measurements. 
These galaxies have moderate SFRs from a few to a 
few tens M yr -1 . The open circles are the galaxies 
with weak near-IR detections, and their measurements 
of SFRs are crude. Because they have large Lya EWs 
and weak UV continuum flux (Figure 6), they show small 
UV SFRs (~ 1 - 2 M yr" 1 ) and high SFR(Lya)-to- 
SFR(UV) ratios. For the purpose of comparison, we 
correct SFR(Lya) for the IGM absorption. Lya emis- 
sion at high redshift is complicated by factors such as 
the resonant scattering of Lya photons and the IGM ab- 
sorption, and it is difficult to model Lya ra diative trans- 
fer and predict intrinsic Lya emission (e.g. , Zheng et al 



2010: lDijkstra et al.llMtllYaiima et al.ll201 2: Silva et al 



20131 ). We thus do not correct individual galaxies. In- 



stead, we derive an average Lya flux reduction for all the 
galaxies. The averag e correction is ca l culate d from the 
composite Lya line of Kashika wa et al.l (|2Qllh by assum- 
ing that the intrinsic Lya profile is symmetric around the 
peak of the composite line. The dotted line in Figure 7 
illustrates the average correction of 24%. When the ef- 
fect of the IGM absorption is taken into account in this 
way, the UV and Lya SFRs for most luminous galaxies 
(filled circles) are reasonably consistent within a factor 
of two. 

Figure 7 shows a trend that the SFR(Lya)-to- 
SFR(UV) ratio becomes increasingly smaller towards 
higher SFRs. The dashed line illustrates this trend 
by displaying the best log-linear fit to the data points. 
This trend is already suggested by the relation between 
EW(Lya) and M1500 in Figure 6, since SFR(UV) re- 
flects M1500 and EW(Lya) mimics SFR(Lya)/SFR(UV). 
Accordingly, this trend has also been shaped by the se- 
lection effect due to the nature of the flux-limited sam- 
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pie, i.e., we have missed low SFR(UV) galaxies with low 
SFR(Lya)-to-SFR(UV) ratios. Therefore, the real trend 
of the SFR(Lya)-to-SFR(UV) ratio, if exists, should be 
much weaker than it appears. This trend could also 
be due partly to the relation between SFRs and the es- 
cape f ractions of Lya photons. The model of lGarel et al.l 
(|2Q12f ) predicts that at high redshift, the escape fractions 
of Lya photons in low-SFR galaxies are close to 1, but 
beco me smaller in galaxies wi th higher SFRs (see also 
e.g., Forero- Romero et al. 2012). This is also consistent 
with I Curt is- Lake et all (|2012h , who found that in a sam- 
ple of luminous (high SFRs) z > 6 LBGs, the Lya SFRs 
are only - 40% of the UV SFRs. 

5. DISCUSSION 
5.1. EW(Lycx) - M 1500 Correlation for LAEs 

Figure 6 shows an apparent correlation between 
EW(Lya) and M1500 in our sample. As we already men- 
tioned, the relation is affected by the nature of the flux- 
limited sample, i.e., the limiting Lya line flux (and there- 
fore the limiting line luminosity) associated with both 
narrow-band and spectroscopic observations. In Fig- 
ure 8(a), we plot the current LAE samples at z ~ 5.7 
and 6.5, as well as those from another lar g e LA E sur- 
vey at z ~ 3.1, 3.7, and 5.7 bv lOuchi etHI (|2008f ). The 
diagonal dotted line is defined by a Lya luminosity of 
2.5 x 10 42 erg s _1 , which roughly corresponds to the lim- 
its of these surveys. The figure illustrates that the slope 
of the EW(Lya)-Muv relation is largely shaped by the 
limiting luminosity, which censors sources that would fall 
below the diagonal dotted lin e. In addition, as pointed 
out bv lCiardullo"eTaII ([20121 ). EW is a derived quantity 
from the emission line and underlying continuum, so any 
relation between EW and the line (or continuum) flux 
may suffer from correlated errors. Such correlated errors 
could result in an appa rent relation by scattering objects 
towards one direction ([Ciardullo et al.|[2012f ). 

Although nothing definite can be said about the prop- 
erties of high-redshift LAEs below our detection limit, 
it is natural to expect that the distribution of LAEs ex- 
tends toward smaller EW(Lya). Therefore, it would be 
interesting to see how the population of low-redshift low- 
luminosity LAEs compares with the high-redshift popu- 
lation shown in Figure 8(a). In Figure 8(b), we include 
the sampl e of GALEX-selected z ~ 0. 2 - 0.4 LAEs pre- 
sented by ICowie. Barger. fe Hul (I20T0I . l20TTbh . The fig- 
ure shows that this low-redshift LAE sample contains 
the kind of LAEs that would be missed at high redshift 
(i.e., below the dotted line). Once these low-luminosity 
LAEs are included, the EW(Lya)-Muv correlation could 
diminish. This is consistent with the fact that no evo- 
lution has been established with t he EW distribution of 
LAE s from z ~ 0.2 - 0.4 to z ~ 3 (ICowie. Barger fc Hul 
[20lQh and from z ~ 3.1 to 5.7 (jOuchi et all 12008( 1. This 
may suggest the possibility that although bright LAEs 
are more abundant at high redshift, the underlying Lya 
EW distribution may be fairly invariant. However, we 
also caution that the lack of a correlation in Figure 8(b) 
is partly caused by another selection effect in these LAE 
samples, namely EW(Lya) > 20 A, which produces a 
sharp horizontal boundary at the bottom. In addition, 
the mix of different samples in Figure 8 could dilute in- 
trinsic relations (if exist). 



Given these selection effects, it is difficult to deter- 
mine with the current LAE samples how strong the cor- 
relation between EW(Lya) and Muv is. To perform 
such an analysis in a statistically meaningful manner, 
we will need a much la rger LAE sample, as pointed out 
by iNilsson et al.l (|2009[ ) . For the discussion in the sub- 
sequent sections, we will assume two extreme cases: (1) 
the slope of the EW(Lya)-Muv correlation is as steep 
as seen in Figure 6 (the dashed line in the lower panel); 
(2) the slope of the EW(Lya)-Muv correlation is com- 
pletely flat. As we shall see, these two limiting cases will 
minimize/maximize the size of the LAE population and 
therefore their contribution to the rest-frame UV lumi- 
nosity density. 

5.2. LAEs and LBGs 

In this paper we call galaxies found by the narrow-band 
technique LAEs and those found by the dropout tech- 
nique LBGs. This is a widely used definition. Strictly 
speaking, this LAE/LBG classification only reflects the 
methodology that we employ to select galaxies. It does 
not necessarily mean that the two types of galaxies are 
intrinsically different. A galaxy with strong Lya emis- 
sion could be detected by the both techniques. Another 
popular definition of LAEs is based on the rest-frame EW 
of the Lya emission line. One galaxy is a LAE if its Lya 
EW is greater than, for example, 20 A. With this defi- 
nition, almost all the galaxies in our sample are LAEs, 
as seen from Figure 6. This definition with Lya EW is 
physically more meaningful. However, the measurements 
of Lya EWs are usually accompanied with large errors. 
Furthermore, one can easily obtain a flux-limited sample, 
but not a EW-limited sample. 

It is not yet totally clear whether high-redshift LAEs 
and LBGs represent physically different populations. Di- 
rect comparison between LAEs and LBGs is difficult. 
The procedure of obtaining a spectroscopic sample of 
LAEs is relatively straightforward. LAE candidates are 
selected based on their Lya luminosities (and one or 
more broad-band photometry), and follow-up spectro- 
scopic identification is also based on their Lya lumi- 
nosities. This results in a complete flux-limited sam- 
ple in terms of Lya luminosity. For LBGs at z > 6, 
candidate selection is based on broad-band colors, but 
spectroscopic identification is based on Lya luminosity. 
Therefore, the resultant LBG sample is inhomogeneous 
in depth of either Lya luminosity or UV continuum lumi- 
nosity, and represents only a subset of LBGs with strong 
Lya emission. 

So far, we did not find significant differences between 
our LAEs and LBGs in the UV luminosity range of 
^1500 < —19.5 mag. Figures 3 and 4 have shown that 
they have similar UV continuum slopes. The LAEs do 
not exhibit steeper slopes than the LBGs, indicating that 
their underlying stellar populations are not very differ- 
ent. To examine the relation between LAEs and LBGs 
further, one useful tool is the EW(Lya)-Muv plot in the 
form of Figures 6 and 8. We expect the two populations 
to have a large overlap, but the focus here is on the non- 
overlapping population(s) that can be picked up by one 
selection method but not by the other. For example, 
the LAE selection may pick up galaxies with extremely 
large EW(Lya) whose continuum emission is so faint that 
they will drop out of LBG samples. Alternatively, the 
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LBG selection can pick up galaxies with small EW(Lya) 
which do not show up in LAE samples. Therefore, we 
would like to know how significant /insignificant such 
non-overlapping populations between LAEs and LBGs 
are. 

Figure 9 plots EW(Lya) and Muv of high-redshift 
LAEs (from Figure 8) and LBGs together. The lat- 
ter sample come s from this work (z ~ 6) and from 
IStark et al.l (j2010D [z ~ 3 - 6). The figure suggests that 
there is no significant difference between the LAE and 
LBG populations in terms of the EW(Lya) and Muv 
distributions. The two populations occupy roughly the 
same region on the plot. For example, we do not see any 
sign of enhanced Lya strengths among the LAE popu- 
lation. This implies that the LBG selection would al- 
most fully recover the LAE population. In other words, 
there are very few LAEs with exceptionally large Lya 
EWs that would drop out of the LBG selection due to 
their faintness in continuum. The reverse (i.e., LBGs 
that would escape the LAE selection due to weak Lya 
emission) is difficult to assess with our current sample 
because our spectroscopic program was not designed as 
an extensive follow-up of LBGs. In the next section, we 
will construct the UV continuum LF of LAEs and ex- 
amine how it compares with the UV continuum LF of 
LBGs. 

5.3. UV Continuum Luminosity Function of LAEs 

Because of the deep HST images, we were able to de- 
tect almost all the LAEs at a significance of > 2<j. This 
allows us to derive the UV continuum LF of LAEs di- 
rectly from the data. The UV LF of LAEs together with 
the UV LF of LBGs will help us estimate the fraction of 
galaxies that have strong Lya emission, and constrain the 
contribution of LAEs to the total UV ionizing photons. 
Our LAEs are from a Lya flux-limited (not EW-limited) 
sample, so they have different detection limits of Lya 
EWs at different UV luminosities (see Figures 8 and 9). 
In this subsection, we will derive the UV LF of LAEs 
with Lya EWs greater than 20 A, by extrapolating the 
observed LAE population down to this EW threshold. 

We first compute the number densities of the LAEs 
without extrapolating the LAE population to the EW 
threshold of 20 A. The area covered by our LAEs is cal- 
culat ed from the area of the whole LAE sample given 
in IKashikaw a et al.l (|2011[ ) , by matching the number of 
our LAEs to the number of the whole LAE sample. We 
do not use the actual area that our HST observations 
covered, because we selected high surface-density regions 
of galaxies for the HST programs (see Section 2.3). 
We then incorpo r ate th e completeness corrections from 
IKashikawa et al.l ()2Q1 lh . Completeness is calculated for 
individual galaxies, as it is a function of narrow-band 
magnitude. The resultant number densities are shown 
as the thick blue and red lines in Figure 10. They set the 
absolute lower limits on the UV LF of LAEs at z = 5.7 
and 6.5. 

We then extrapolate the densities to EW = 20 A using 
the Lya EW distribution function. Figure 11 shows the 
Lya EW distribution of the LAEs in our sample. We 
assume that the EW distribution is an exponential func- 
tion n oc e - W w °, where w is EW and wq is the e-folding 
width. It is often assumed that wq is correlated with 



UV l uminosity (e.g. Kashika wa et al.l 120111 : IStark et al.l 
l20lTh . We assume that wq is a log-linear function of 
^1500 : log(^o) = a + bx (M1500 + 2O), where a and b are 
scaling factors. Given the small numbers of the galax- 
ies, it is not realistic to obtain reliable measurements for 
both a and b. Therefore, we consider the two extreme 
cases mentioned in Section 5.1: (1) b = 0.225, where we 
assume that the slope of the EW(Lya)-Mi5oo correla- 
tion is as steep as seen in Figure 6 (the dashed line in 
the lower panel); (2) 6 = 0, where the slope of the cor- 
relation is completely flat. These two limiting cases will 
minimize (b = 0.225) /maximize (b = 0) the size of the 
LAE population. The real b value is between and 0.225. 
We then determine a by fitting the exponential function 
to the EW distribution of the LAEs with Mi 500 < -19.5 
mag in our sample. The best fit is shown in the top panel 
of Figure 11. With the derived a and assumed b we can 
calculate the Lya EW distribution at any luminosity. 
The other panels in Figure 11 show the predicted EW 
distributions for the two extreme cases with the dashed 
(b = 0.225) and dash-dotted (b = 0) lines. 

The shaded regions in Figure 10 show the UV LFs of 
LAEs corrected to EW(Lya) = 20 A for the two extreme 
cases. Their lower and upper boundaries represent the 
lower and upper limits of the LFs. The la statistical 
uncertainties have been included. At M1500 < —20.5 
mag, no correction is actually applied, because our sam- 
ple is complete at EW = 20 A down to this magnitude. 
Towards fainter magnitudes, the correction becomes in- 
creasingly larger, and the difference of correction between 
the two cases also becomes larger. In the faintest end, 
the difference of the number densities between the lower 
and upper limits is larger than a factor of 3. 

The cyan and green curves in Figure 10 show 
the U V LFs of photo-sele c ted LBGs from the liter- 
ature dBouwens et al.l I2007L 120111 : iMcLure et all 120091 : 
lOuchi et al.l l2QQ9h . compared to our results. At the 
bright end, the LAE UV LFs are roughly comparable to 
th e LBG LFj^l cons isten t with similar results re ported 
bv lOuchi et "all (|2008[ ) and IKashikawa et al.l (|2011f ). This 
suggests a large fraction of LAEs among the bright- 
est LBGs. Such high fraction has been observed by 
iCurtis-Lake et al.l (|2012l ). who found that the LAE frac- 
tion in a sample of very bright (Muv < —21 mag) UDS 
LBGs is about 50%. IStark et al.l (j2011f ) reported a lower 
fraction (20 ± 8%) of LAEs with EW(Lya) > 25 A in a 
sample of bright z~6 LBGs (-21.75 < M uv < -20.25). 
We note, however, their sample contains a small fraction 
(10%) of LAEs with M uv < -21 mag, so their LAE 
fraction is dominated by the galaxies at Muv > — 21 
mag. Our LAE fraction (~ 35%) in a similar range of 
—21 < Muv < —20 is not much different from theirs, 
given their slightly higher EW limit. Note that the 
measurements of the LAE fraction among the brightest 
galaxies are subject to large uncertainties owing to the 
small numbers of galaxies available. 

In Figure 10, the UV LF of LAEs at the faint end is 
even more uncertain, as it depends critically on how we 
extend the EW(Lya) distribution toward EW = 20 A, 

11 Strictly speaking, our first HST GO program 11149 slightly 
favored galaxies with bright continuum flux. This could increase 
the number densities by up to 50% in the bright end, but is not 
large enough to change the basic conclusion. 
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which is far below our detection limits. Current ground- 
based spectroscopy is not yet able to reach this regime. 
Depending on the value of b we assume, the faint-end 
slope of the LAE UV LF could be almost as steep as 
that of the LBG UV LF (b - 0) or significantly flatter 
(b > 0). If we maximize the LAE UV LF by assuming 
b ~ 0, we have to apply a large incompleteness correction 
(> xlO at the faintest end) to account for LAEs down 
to EW(Lya) = 20 A. It means that the Lya EW dis- 
tribution must be highly concentrated on small EWs, so 
that the observed LAEs represent a small sub-group of 
high-EW galaxies that constitutes a tiny fraction of the 
overall LAE population. Despite this uncertainty, it is 
reasonable to conclude with Figure 10 that the number 
density of LAEs falls significantly short of that of LBGs 
at the faint end. This result could have an important 
implication for the population of galaxies that were re- 
sponsible for cosmic reionization. It implies the existence 
of a large population of LBGs with weak Lya emission 
(EW < 20 A). Because the UV LF slope of LBGs is steep, 
the vast majority of the UV photons would come from 
very faint galaxies (far below L*). This suggests that the 
LBG population with weak Lya emission dominate the 
UV photon budget for cosmic reionization. 

In Figure 10, the comparison between the UV LFs of 
LAEs at z ~ 5.7 and 6.5 is straightforward, and it shows 
little evolution. This supports the previous claim that 
the Lya LF evolution from z ~ 5.7 to 6.5 is due to an 
increasing neutral fraction of the IGM between the two 
redshifts an d not due to the evolut i on of the LAE galaxy 
population (Kashika wa et al.ll2006l l201lt ). These studies 
showed that the Lya LF evolves rapidly from z ~ 5.7 to 
6.5, and pointed out in particular that there is a lack of 
luminous LAEs at z ~ 6.5. In contrast, the evolution of 
the UV LF of LAEs between the two redshifts is much 
more modest. This was explained by the increasing neu- 
tral fraction of IGM from z ~ 5.7 to 6.5 that attenuates 
Lya emission. In these studies, the UV continuum lu- 
minosities were measured from the SDF z'-band imaging 
data, which can only detect bright LAEs, and are also 
contaminated by Lya emission (and Lyman break) for 
z ~ 6.5 LAEs. Our HST data were deep enough to de- 
tect almost all the LAEs in our sample, so our measure- 
ments of UV continuum luminosities are more robust. 
We confirm the consistency of the UV LFs of LAEs at 
z ~ 5.7 and 6.5, which rules out the possibility of a strong 
intrinsic galaxy evolution between the two redshifts, and 
strengthens the interpretation that the increasing neu- 
tral fraction of IGM causes the strong evolution of the 
Lya LF towards z ~ 6.5. 

6. SUMMARY 

We have carried out deep HST near-IR and Spitzer 
mid-IR observations of a large sample of spec-confirmed 
galaxies at 5.6 < z < 7. The sample contains 51 LAEs 
and 16 LBGs, representing the most luminous galaxies in 
this redshift range. The LAEs were from a complete flux- 
limited sample. The LBGs have quite different depth, 
and only contain those with strong Lya emission. The 
majority of the galaxies (62 out of 67) were discovered in 
the SDF, and the remaining were found in the SXDS. The 
HST observations of the SDF galaxies were made with 
WFC3 and NICMOS in the J (F110W or F125W) and H 
(F160W) bands. The depth is two HST orbits per band 



for most of the objects. With such depth nearly 80% of 
the galaxies were detected at high significance (> 5<r) in 
the J band. The Spitzer observations of the SDF were 
made in two IRAC channels 1 and 2. The depth varies 
across the field from ~3 hrs to ~6 hrs. The infrared 
data of the five SXDS galaxies were taken from the HST 
and Spitzer archive. In addition to the infrared data, 
we also have extremely deep optical images in a series 
of broad bands (BVRi'z'y) and narrow bands (NB816, 
NB921, and NB973). We have used the combination of 
the optical and infrared data to derive the properties of 
rest-frame UV continuum and Lya emission, such as UV 
continuum luminosities and slopes, Lya luminosities and 
EWs, and SFRs etc. 

While the whole sample covers a large UV continuum 
luminosity range from M1500 < —22 mag to M1500 > — 18 
mag, the galaxies with significant detections (> 5a) in 
the J band cover a bright range of M1500 < —19.5 mag. 
These galaxies have steep UV continuum slopes roughly 
between j3 c± —1.5 and j3 —3.5, with a weighted mean 
of /3 c± —2.3. This value is slightly steeper than the slopes 
of photo-selected LBGs reported in previous studies in 
the literature. This is due to the fact that our galaxies 
have strong Lya emission, and thus lower dust extinction 
and steeper UV slopes. The slope j3 shows little trend 
with the UV luminosity when the several brightest galax- 
ies are excluded. The LAEs do not display significantly 
bluer slopes than the LBGs in this sample and those in 
previous studies, suggesting that most LAEs are proba- 
bly not younger than LBGs as expected previously. A 
small fraction of our galaxies have extremely steep UV 
slopes with j3 c± — 3. Stellar populations in these galax- 
ies could be very young with extremely low metallicity 
and dust content. Our galaxies have moderately strong 
Lya EWs in a wide range of -10 to -200 A. The SFRs 
estimated from the Lya and UV luminosities are also 
moderate, from a few to a few tens M yr _1 . 

We have also derived the UV LFs of LAEs with EW > 
20 A at z ~ 5.7 and 6.5, using our deep HST near- 
IR data. We have confirmed that the number density 
of LAEs (EW > 20 A) is comparable to that of LBGs 
at the bright end. This suggests that the LAE fraction 
among bright LBGs is very high. We also concluded that 
the number density of LAEs cannot be as high as that of 
LBGs at the faint end, even if we maximize the former by 
adjusting a model parameter b (i.e., 6 = 0). This implies 
that there exists a substantial population of faint LBGs 
with weak Lya emission (EW < 20 A) that could be 
the dominant contribution to the total UV luminosity 
density at z > 6. 

The LAEs and LBGs in this sample are indistinguish- 
able in many aspects of the Lya and UV continuum prop- 
erties. In fact, in the EW(Lya)-Muv diagram, the dis- 
tributions of the two populations are quite similar, show- 
ing no sign of enhanced Lya strengths among the LAE 
population. All these seem to indicate that LAEs can be 
considered as a subset of LBGs with strong Lya emission 
lines. 
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APPENDIX 

THUMBNAIL IMAGES OF THE GALAXIES 

Thumbnail images of the galaxies in one of the three Subaru Suprime-Cam narrow bands (NB816, NB921, and 
NB973, depending on redshift), three Suprime-Cam broad bands iz'y, two HST near-IR bands J125 (or Jno) and 
-H16O5 an d two Spitzer IRAC bands. The galaxies are in the middle of the images. The size of the images is 6 / /6 x 6 / /6 
(north is up and east to the left). 
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Note. — The first 62 galaxies are from the SDF, and the last 5 galaxies are from the SXDS. 

a NB = NB816 fo r z ~ 5.7 LAEs, NB = N B921 for z ~ 6.5 LAEs , NB = NB9 73 for the z ~ 7 LAE. 

b Ref erenc es: (1) [Shimasaku et al.l (f2006h; (2) IKashikawa et al.l (f20lTl); ( 3 ) IQta et al.l (f2Q08h; (4) INagao et al.l (120071): (5) INagao et al.l 
(2005); (6) Jiang et al. (2011); (7) Nagao et al. (2004); (8) Taniguchi et al. (2005); (9) Kashikawa et al. (2006); (10) lye et al. (2006); (11) 
IKodaira et al. (2003); (12) Curtis-Lake et al. (2012); (13)[(7uchi et al. (201Cj); (14)[Ouchi et al. (200^); (*) unpublished. 
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TABLE 2 

Spectral Properties of the Galaxies 



No. 


M1500 


UV Slope 


L(Lya) 


EW(Lya) 


SFR(UV) 


SFR(Lya) 






(0) 


CIO 42 erg s _1 ) 


(A) 


(Mm yr _1 ) 

V ±v - L W J 1 - ) 


(Mm yr _1 ) 

V ±V - L W J 1 - ) 


2 


-19.53 ±0.15 


-2.09 ±0.52 


13.5 ± 1.2 


231.3 ±38.6 


3.8 ±0.5 


12.3 ± 1.1 


3 


-21.46 ±0.03 


-2.36 ±0.13 


10.5 ±0.5 


28.6 ± 1.6 


20.9 ±0.5 


9.5 ±0.5 


4 


-20.50 ±0.27 


-1.96 ±0.79 


8.7 ±0.5 


62.7 ± 16.0 


9.9 ±2.5 


8.0 ±0.4 


5 


-20.07 ±0.09 


-3.00 ±0.35 


8.7 ±0.6 


75.2 ±8.0 


4.9 ±0.4 


7.9 ±0.5 


7 


-20.12 ±0.08 


-3.01 ±0.34 


7.3 ±0.5 


60.2 ±5.9 


5.1 ±0.4 


6.7 ±0.4 


10 


-19.69 ±0.13 


-2.89 ±0.49 


2.7 ±0.4 


33.8 ±6.3 


3.5 ±0.4 


2.5 ±0.4 


14 


-21.27 ±0.04 


-2.19 ±0.28 


6.4 ±0.5 


21.6 ± 1.9 


18.4 ±0.7 


5.8 ±0.5 


15 


-21.94 ±0.02 


-2.21 ±0.11 


20.8 ±0.8 


37.8 ± 1.7 


33.9 ±0.6 


19.0 ±0.7 


17 


-19.99 ±0.34 


-3.02 ±1.10 


3.0 ±0.4 


28.0 ±9.7 


4.5 ± 1.4 


2.8 ±0.4 


19 


-20.46 ±0.08 


-1.57 ±0.30 


4.0 ±0.5 


32.0 ±4.4 


10.9 ±0.8 


3.6 ±0.4 


20 


-20.46 ±0.08 


-1.68 ±0.26 


16.5 ±0.9 


131.1 ± 11.4 


10.5 ±0.7 


15.1 ±0.8 


21 


-19.85 ±0.10 


-3.18 ±0.37 


4.4 ±0.6 


44.9 ± 7.6 


3.8 ±0.3 


4.0 ±0.6 


22 


-19.90 ±0.12 


-2.24 ±0.44 


5.1 ±0.7 


60.9 ±10.4 


5.1 ±0.6 


4.7 ±0.6 


23 


-21.17 ±0.03 


-2.35 ±0.14 


5.0 ±0.8 


18.0 ±2.8 


16.0 ±0.5 


4.6 ±0.7 


24 


-21.19 ±0.04 


-2.58 ±0.22 


13.6 ±0.8 


45.6 ±3.1 


15.2 ±0.5 


12.4 ±0.7 


25 


-21.30 ±0.04 


-2.82 ±0.20 


14.6 ± 1.2 


42.0 ±3.8 


15.9 ±0.6 


13.3 ± 1.1 


26 






17.4 ± 1.2 






15.8 ± 1.1 


27 


-20.20 ±0.09 


-2.20 ±0.40 


6.9 ±0.8 


59.5 ±8.6 


6.9 ±0.6 


6.3 ±0.7 


28 


-20.85 ±0.06 


-2.81 ±0.39 


18.3 ± 1.2 


80.7 ±7.1 


10.4 ±0.6 


16.6 ± 1.1 


29 


-19.58 ±0.12 


-2.62 ±0.47 


2.9 ±0.8 


41.3 ±12.6 


3.4 ±0.4 


2.6 ±0.7 


30 


-20.68 ±0.07 


-3.45 ±0.38 


4.6 ±0.8 


20.5 ±3.9 


7.8 ±0.5 


4.2 ±0.8 


31 


-20.53 ±0.19 


-2.51 ±0.76 


18.8 ±0.4 


117.4 ±20.2 


8.4 ± 1.4 


17.1 ±0.4 


32 






5.9 ±0.8 






5.3 ±0.8 


33 


-20.03 ±0.24 


-2.40 ±0.95 


4.4 ±0.9 


44.7 ±13.4 


5.5 ± 1.2 


4.0 ±0.8 


34 


-21.82 ±0.04 


-1.79 ±0.17 


4.1 ±0.9 


8.6 ±1.9 


35.0 ±1.4 


3.7 ±0.8 


35 


-21.37 ±0.11 


-2.31 ±0.41 


18.1 ±0.5 


54.2 ±5.5 


19.5 ±1.9 


16.5 ±0.4 


36 


-21.25 ±0.07 


-2.57 ±0.30 


2.8 ±1.0 


8.9 ±3.1 


16.1 ± 1.0 


2.5 ±0.9 


37 


-20.19 ±0.26 


-2.55 ±0.94 


12.6 ± 1.0 


106.5 ± 26.6 


6.1 ± 1.4 


11.5 ±0.9 


39 


-20.21 ±0.20 


-1.69 ±0.88 


16.1 ±0.8 


161.1 ±31.3 


8.3 ± 1.6 


14.7 ±0.8 


40 


-20.39 ±0.12 


-4.37 ±0.73 


1.4 ±0.2 


6.8 ±1.2 


5.2 ±0.6 


1.3 ±0.2 


43 


-20.70 ±0.09 


-3.05 ±0.47 


5.2 ±0.4 


24.6 ±2.9 


8.6 ±0.7 


4.7 ±0.4 


44 


-21.74 ±0.03 


-1.79 ±0.16 


19.9 ±0.8 


47.4 ±2.3 


32.6 ±1.0 


18.1 ±0.7 


45 


-19.99 ±0.31 


-3.33 ±1.57 


3.4 ±0.4 


28.9 ±9.2 


4.2 ± 1.2 


3.1 ±0.4 


46 


-20.13 ±0.21 


-1.87 ±0.91 


4.1 ±0.5 


42.7 ±9.5 


7.2 ± 1.4 


3.8 ±0.4 


47 


-22.08 ±0.03 


-1.57 ±0.17 


5.4 ±0.5 


9.8 ±0.9 


48.4 ± 1.5 


4.9 ±0.4 


48 


-19.64 ±0.30 


-2.09 ±1.39 


5.2 ±0.5 


80.8 ±24.2 


4.3 ± 1.2 


4.8 ±0.5 


49 


-20.49 ±0.13 


-1.98 ±0.59 


5.9 ±0.5 


42.5 ±6.3 


9.7 ± 1.1 


5.4 ±0.5 


50 


-20.09 ±0.19 


-1.60 ±0.82 


9.6 ±0.6 


108.7 ±20.2 


7.7 ± 1.3 


8.7 ±0.6 


52 


-19.31 ±0.24 


-2.44 ± 1.14 


6.3 ±0.6 


121.6 ±29.1 


2.8 ±0.6 


5.7 ±0.5 


54 


-21.12 ±0.09 


-2.64 ±0.42 


13.0 ±0.7 


45.6 ±4.3 


14.1 ± 1.1 


11.8 ±0.6 


58 


-21.08 ±0.09 


-1.84 ±0.40 


5.3 ±0.5 


23.1 ± 2.8 


17.5 ±1.4 


4.9 ±0.4 


61 


-21.06 ±0.08 


-3.47 ±0.48 


12.6 ±0.7 


39.3 ±3.5 


11.0 ±0.8 


11.4 ±0.6 


62 


-21.51 ±0.06 


-2.09 ±0.31 


15.1 ± 1.6 


41.5 ±5.0 


23.9 ±1.4 


13.7 ± 1.4 


63 


-21.67 ±0.04 


-2.78 ±0.17 


18.2 ±3.0 


37.7 ±6.3 


22.4 ±0.8 


16.5 ±2.7 


64 


-21.62 ±0.10 


-3.01 ±0.38 


13.0 ± 1.1 


26.6 ±3.3 


20.2 ±1.9 


11.8 ± 1.0 


66 


-20.10 ±0.31 


-1.82 ±1.35 


8.3 ±1.0 


89.2 ±27.5 


7.1 ±2.0 


7.6 ±0.9 


67 


-22.08 ±0.08 


-1.87 ±0.33 


22.9 ± 1.2 


39.3 ±3.5 


43.3 ±3.1 


20.8 ± 1.1 
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0.6 0.8 1.0 1.2 1.4 1.6 
Wavelength (/xm) 

Fig. 1. Transmission curves of six Subaru Suprime-Cam filters {%' ', z', y, NB816, NB921, and NB973) and three HST WFC3 filters 
(F110W, F125W, and F160W). System responses such as detector quantum efficiencies are included. Note three narrow-band niters NB816, 
NB921, and NB973 (dotted black curves) at 0.8 ~ 1.0 /mi. 
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Fig. 2. — Examples of measuring UV continuum and Lya line emission. The upper panel illustrates our model fit to a LAE at z = 5.710. 
The red circles are the photometric points that we use. The broad-band data in the z'yJuoHiQo bands are used to fit the power-law 
continuum (see Eq. 1). The horizontal bars indicate the wavelength coverage of the broad bands. The photometric error bars are also 
shown, but the errors are very small. When IGM absorption is applied, the continuum flux at the blue side of the Lyo: line is absorbed 
(the dashed line). The region around the Lyo: emission line (the dotted line) is zoomed in to show this absorption. Finally the factor A in 
Eq. 1 is calculated by scaling the Lyo; emission line to match the NB816-band photometry. The solid profile is the final model spectrum 
(continuum + Lyo;) for this galaxy. The lower panel shows our model fit to a LAE at z = 6.543. In this case we use the data in the 
2/^125^160 bands to compute its continuum, and use the NB921-band photometry to measure its Lyo; emission. 
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-22 -21 -20 -19 -5 -4 -3 -2 -1 

M 1500 Slope (jS) 

Fig. 3. — (a) UV continuum slope f3 as a function of the UV continuum luminosity Misoo- The blue and red circles represent the LAEs 
at z ~ 5.7 and 6.5 (including z ~ 7), respectively, and the green circles represent the LBGs at z ~ 6. The slopes in this sample have a 
weighted mean (3 ~ —2.3. The dashed line is the best linear fit to all the data points. The dash-dotted line is the best fit to the galaxies 
fainter than M1500 = —21.7 mag, which suggests that the slope (3 shows little trend with M1500 when the several brightest galaxies at 
M1500 ~ — 22 mag are excluded, (b) Normalized distribution of the observed /3 (histogram) compared with the distribution of (3 expected 
(the red profile) if all the galaxies have an intrinsic (3 = —2.3 (the vertical dotted line) with the observed uncertainties. The expected (3 
distribution is the average of the Gaussian distributions for individual galaxies. The observed distribution is broader and flatter, meaning 
significant intrinsic scatter of (3. It clearly shows a statistically significant excess of galaxies with f3 ~ — 3, indicating that the existence of 
f3 ~ — 3 galaxies are statistically robust. 
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Fig. 4. — UV continuum slope /3 as a function of the UV continuum luminosity M1500, compared to previous studies. The gray points with 
error bars represent the galaxies shown in Figure 3(a). The magenta star represents the stacked image of the 9 objects with crp > 0.8 that 
have the J125 and ifi60-band images. The magenta square represents the stacked image of 10 galaxies with f3 < —2.8 that have the J125 
and .ffi60-band images. The horizontal bars on the magenta symbols indicate their magnitude ranges. The blue crosses, g reen trian gles, 
and red diamonds display the average slopes from previous studies: the blu e crosses are the average (3 at z = 5 — 7 from Dunlop et al. 
(2012a); t he gree n triangles are the average (3 at z ~ 6 from Bouwens et al. (2012b); the red diamonds are the average (3 at z ~ 6 from 
Finkelstein et al. (2012). These points are connected by solid lines to guide eyes. At M1500 < —19.5 mag, these studies do not show 
evidence of a significant correlation between M1500 and (3. Such relation, if exists, could be driven by very faint galaxies. 
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Fig. 5. — Comparison between the observed Lya luminosities of the LAEs in our paper and the Lya luminosities given in Kashikawa et al. 
(2011). The upper panel compares our luminosities to those derived from the photometric data in the galaxy discovery papers. The 
photometric data used in these papers are usually the narrow-band and z'-band photometry. The blue and red circles represent the LAEs 
at z ~ 5.7 and 6.5, respectively. The dotted line is the locus where the two measurements are equal. The error bars are typical measurement 
errors. The two measurements agr ee well with ea ch other. The lower panel shows the comparison between our results and the results from 
the spectroscopic data in Kashikawa et al. (2011). The two are also consistent, with a slightly larger scatter than that in the upper panel. 
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Fig. 6. — Lyo; luminosity and rest-frame EW (A) as a function of UV luminosity Misoo- The color-coded filled circles have the same 
meaning as in Figure 3. The open circles represent the LAEs with weak detections (< 5cr) in the J band, and their measurements of 
physical quantities are crude (see Section 4.2). The dashed lines are the best linear fit to the filled circles. The Lyo; luminosity shows weak 
dependence on M1500 in the upper panel. The lower panel shows that lower-luminosity galaxies tend to have higher Lyo; EWs. The real 
relation should be weaker than it appears due to the nature of the flux-limited sample. There is no significant difference between LAEs 
and LBGs in the range of M1500 < —19.5 mag. 
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Fig. 7. — SFRs estimated from UV continuum luminosities and Lya: luminosities in units of Mq yr~ 1 . The color-coded symbols have 
the same meaning as in Figure 6. The solid line represents the position where the two measurements are equal. The dotted line is the 
expected ratio of the Lya SFRs to the UV SFRs assuming that Lya is absorbed by 24% on average. The dashed line is the best linear 
fit (in logarithmic space) to the relatively robust measurements (filled circles). The Lya SFRs are increasingly smaller than the UV SFRs 
towards higher SFRs, mostly due to the selection effect shown in Figure 6. It could also be due partly to the smaller escape fractions of 
Lya photons in higher-SFR galaxies. 



1000 F 



(a) LAEs: High-z 



100: 



< 



10 r 



•ft 



• LAE z=5.7 (this work) 

• LAE z=6.5 (this work) 
LAE z=3.1 (Ouchi08) 

• LAE z=3.7 (Ouchi08) 

• LAE z=5.7 (Ouchi08) 



V 



-24 



-23 



-22 



-21 



-20 



-19 



-18 



-17 



1000 F 
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Fig. 8. — EW(Lyo;) vs. Mjjy for LAE s. (a) EW(Lya) is p lotted as a function of Mtjv for the z ~ 5.7 and 6.5 LAEs in our sample, as well 
as for the z ~ 3.1, 3.7 , and 5.7 LAEs by Ouchi et al. (2008). As for Muv, we used M1500 listed in Table 2 for our sample. For the sample 
of [Quchi et al. (2008), we used the published Muv, which correspond to M1300-I6OO depending on redshift. The diagonal dotted line is 
defined by a limiting Lya line luminosity of 2.5 x 10 42 erg s _1 , which roughly corresponds to the limits of these LAE s urveys. The slope of 
the E W(LycQ-Mxjv correlation is largely shaped by the limiting luminosity, (b) GALEX-selecte d z ^ 0.2 — 0.4 LAE s ( Cowie, Barger, & Hu 
2011b) are added for comparison. Muv was derived from the GALEX near-UV photometry in Cowie, Barger, & Hu (2010). The effective 
wavelength of the near-UV band is 2315.7 A (Morrissev et al. 2007), so at redshifts z ~ 0.2 — 0.4, the near-UV band is sampling ~1650-1930 
A in the rest-frame. The near-UV, instead of the far-UV band, was used because the latter contains the redshifted Lya emission. The 
figure indicates that the EW(Lya)-Muv relation is largely shaped by the limiting luminosity. 



23 



1000 F 



60 

a 

< 



High-z LAEs vs. LBGs 



100 




10 

\<* ■ . - 

• LAE z=5.7 (this work) 

• LAE z=6.5 (this work) 
LAEz=3.1 (Ouchi08) - 

• LAE z=3.7 (Ouchi08) 

• LAE z=5.7 (Ouchi08) 
LBGz~6 (this work) 

■ LBG z=3-6 (Stark 10) 

1 [ 

-24 -23 -22 -21 -20 -19 -18 -17 



Fig. 9. — Same as Figure 8, but for LBGs from this work (2 ~ 6) and from Stark et al. (2010) (2 ~ 3 — 6). A number of LBGs from 
Stark et al. (2010) fall below the L(Lyoi) = 2.5 x 10 42 erg s — 1 line, reflecting the deeper limiting Lyo; line flux of the data. The overall 
distribution of the points appears similar between the LAE and LBG samples. 
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Fig. 10. — UV continuum LFs of LAEs with EW > 20 A at z ~ 5.7 and 6.5. The thick blue and red lines represent the spatial densities 
of the LAEs in our sample, without incompleteness correction to EW = 20 A. The shaded regions show the UV LFs of LAEs corrected to 
EW(Lya)= 20 A. Their lower and upper boundaries represent the lower and upper limits of the LFs. The la statistical uncertainties have 
been inclu ded. The cya n and green curves represent the UV LFs of LBGs at z ~ 5.9 and 6.8 JB ouwens et al. 2007, 2011; McLure et al. 
2009; Ouchi et al. 2009). Compared to LBGs, LAEs have comparable spatial densities in the bright end, but significantly lower densities 
in the faint end. 
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Fig. 11. — Normalized Lyo: EW distribution of the LAEs in our sample. The vertical dotted lines represent the approximate EW limits 
for complete LAE samples at the given M1500 ranges. We use an exponential function to fit the EW distribution. We adopt two extreme 
e-folding widths (b = 0.225 and 0) in Section 5.3. The solid curve in the top panel is the best fit to the LAEs at M1500 < —19.5 mag (it is 
the same for the two cases). The dashed (b = 0.225) and dash-dotted (b = 0) curves in the other panels are the predicted EW distributions 
for the two extreme cases. The real distribution should be between the two. 
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Fig. 12. — Thumbnail images of the galaxies in our sample. The galaxies are in the middle of the images. The size of the images is 
6'/ 6 x 6'/ 6 (north is up and east to the left). An image is blank if the data is either not available or meaningless in this band. 
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Fig. 12. — Continued. 
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Fig. 12. — Continued. 



